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ABSTRACT 

We study the Abelian Higgs Model using an improved form of the action in 
the scalar sector. The subleading corrections are carefully analysed and the con- 
nection between lattice and continuous parameters is worked out. The simulation 
shows a remarkable improvement of the numerical performance. 



1 Introduction 



The Abelian Higgs Model has mostly been studied in recent years as a theoretical 
laboratory in the context of the Electroweak Baryogenesis scenario. As it is well 
known by now, the lattice investigations of the model are very demanding in 
computer power and time. It would be helpful to use an improved form of the 
lattice action, to reduce the autocorrelation and come closer to the continuum 
limit. 

We did not try an improvement on the pure gauge sector of the action. The 
strategy which is usually followed does not seem very reliable: actually it leaves 
behind several subleading contributions, so it is not very efficient in eliminating all 
of the unwanted terms. On the other hand, the weak coupling regime where these 
simulations are performed suggest that the improvement of the gauge sector may 
not be very important. We mainly concentrated on improving the scalar sector; 
the procedure suffers more or less from the same problems, however even a modest 
improvement is very important in this case. The model on which we will work 
has been treated before ([|l|]), so we have reference results to compare with. 

A promising approach to study the four- dimensional model at finite tempera- 
ture is through reduction to an effective model in three dimensions. This can be 
done if the couplings are small and the temperature is much larger than any other 
mass scale in the theory. The parameters of the reduced theory are related to the 
ones of the original model through perturbation theory. The reduced theory has 
some advantages over the original one from the computational point of view. It is 
super-renormalizable and yields transparent relations between the (dimensionful) 
continuous parameters and the lattice ones. Moreover, the number of mass scales 
is drastically reduced: (a) the scale T, present in four dimensions is evidently 
absent, (b) one may also integrate out the temporal component Aq of the gauge 
field, so its mass scale gT also disappears. Thus there are two mass scales less 
and this reduces substantially the computer time needed to get reliable results. 

The action for the U(l) gauge-Higgs model at finite temperature is: 

/■/3 /• 1 

S[A^{t,x),^{t,x)]= dr J d^x[-F^^F^^ + \D^ip\^ + m^ifi*ifi + X{ip*ip)% (1) 
where /3 = 1/T. 

If the action is expressed in terms of Fourier components, the mass terms are 
of the type: 

[{2nnTy + ikY]\A,in,k)\' (2) 
[i27rnTr + ikrMn,k)\' (3) 

where n = — oo, . . . , oo. 

At high temperatures T and energy scales less than 27iT the non-static modes 
Afj.{n 0, k), ip{n 0, k) are then supressed by the factor (27rnT)^ relative to the 
static Af^{n = 0, k) and ip{n = 0, k) modes. The method of dimensional reduction 
consists in integrating out the non-static modes in the action and deriving an 
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effective action (0). We notice that the mass of the adjoint Higgs field is of 
order gT, which is large compared to g'^T, the typical scale of the theory. Thus, 
we integrate it out using perturbation theory (@). 

The effective action may then be written in the form: 

S3D eff[Aiix), ^3{x)] = J d^xl^FijFij + I A'/'sr + mly^l^s + Xsiy^ly^sf] (4) 

The index 3 in (|) is for the 3D character of the theory. The relations between 
4D and 3D parameters are (up to 1-loop): 

9l = (5) 



A3 = TA(/.) - (6) 

< = IqHi^W'. (7) 

It is convenient to use the new set of parameters [g"^, x, y) rather than the set 
i.9lAi.^l4)- are defined as (0]): 

x = ^ (8) 

93 



^1(9: 

-,4 
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y = (9) 
93. 

It is evident that x is just proportional to the ratio of the squares of the scalar 
and vector masses; y is related to the temperature. The parameters x, y can be 
expressed in terms of the four-dimensional parameters as follows([^]): 

x = l4^-^ (10) 

We have concentrated on the phase transition line. We have chosen to fix the 
Higgs mass to a fixed value (30 GeV), g to 1/3, m^y to 80.6 GeV and study the 
characteristics of the phase transition. 
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2 The improved lattice action 



The whole idea of improving the lattice action has been put forward (0, [0, |Q, 
[1]) 1131) enhance performance of the lattice calculations. The lattice actions, 
when expanded in powers of the lattice spacing a, yield the terms of the continuum 
action plus subleading terms, i.e. terms multiplied by higher powers of a. The 
procedure is to include additional terms in the action, so that the corrections 
that remain in the naive continuum limit start at a higher power of the lattice 
spacing, as compared to the usual action. This work follows most closely |]10 . 



2.1 The pure gauge part 

The pure gauge part is expressed by the plaquette term in the action. We use the 
non-compact formulation; the initial action is defined by I3g Z]o<i<i -^4'; "^^ere 
Fij = A{Aj{x) - A^jAi{x), A{Aj{x) = Aj{x + i) - Aj{x). In the following we 
treat the part of the action having to do with the xy plane, i.e. we consider a 
two dimensional model; generalization to include the remaining hyperplanes is 
straightforward. In addition, we consider for both the gauge and the scalar field 
part of the action two versions of the improvement: the "continuum" and the 
lattice version. In the former lattice with spacing equal to a is embedded 

in the continuum space-time and objects resembling the usual lattice quantities 
are considered. An interpolation is used, which makes it easy to eliminate all 
of the subleading terms both in the gauge and the scalar field sectors. The 
lattice approach is the treatment of the actual lattice model. It is not in general 
possible to eliminate all of the terms and one should be content with a partial 
cancellation. 



2.1.1 Lattice embedded in the continuum 

We consider the quantities: 



1 I 1 



Cu = J ^ dtA^{x + at,y - -) + J ^ dtAy{x + -,y + at) 

^ dtA,j,{x + at,y + -) - j ^ dtAy{x - -,y + at), 
dtAx{x + at,y — -) + / dtAy{x + a,y + at) 
— / dtAx{x + at,y + a) — / dtAy(x — a, y + at) , 
C21 = J ^ dtAx{x + at,y — a) + J dtAy{x + -,y + at) 
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— dtAx(x + at,y + a) — / dtAy(x — -,y + at). 
J-\ J-i 2 

Notice that Cu represents a continuum version of the 1x1 plaquette, while the 
two terms Cu, C21 represent the 1x2 and 2x1 plaquettes. 

One then expands these quantities in powers of a and ends up with: 

a ^Cii ~ [dxAy — dyAx] + —(i^[dxxxAy — dyyyA^ + dxyyAy — dxxyAx] (12) 

a ^Ci2 — 2[dxAy — dyAx] + -:^(i^\4dxxxAy — dyyyAx + dxyyAy — 4:dxxyAx] (13) 

a'~^C2i ~ 2[dxAy - dyAx] + ^a^idxxxAy - ^dyyyAx + ^dxyyAy - dxxyAx] (14) 

The terms that will appear in the action can be written more simply by using 
the notations: F = dxAy dyAx-, P = dxxxAy dyyyAx + dxyyAy dxxyAx- They 
read, up to O(a^) : 

a-2(Cn)' ^F^ + ^a'FP, a-\C^2)' + a-^C^.f 8F' + L^FP (15) 
Now it is easy to write down the expression for the action up to O(a^) : 

a-'S = Y^[Aa-'\Ci,f + Ba-'{{Cuf + (C21)'}] 
x,y 

~ Y.[A{F^ + ^a^FP) + B{8F^ + ^a^FP)]. (16) 

x,y 1^ 3 

We observe that in general we get the continuum action plus the FP terms, 
which are lattice artifacts. If our aim is to better approach the continuum action, 
we should arrange that these artificial terms vanish; in addition, the coefficient 
of the term should be one. This is easy in this approach: we just choose 
^4 = 1, B = Thus the improved action reads: 

E [l{C^,'{x,y)y-^iC^^ix,y)r-l-iC^i:i^,y)r] 

x,y,iJ.<u 1^ 1^ 



2.1.2 Actual lattice formulation 

In actual lattice calculations one cannot use the interpolation of the previous 
section, namely the one based on the t-integrations. One has link variables and 
the plaquettes used in the action are sums of four (or six) such variables. It is 
therefore very interesting to find out how the above analysis is modified if the 
real situation on the lattice is considered. We will use the same plaquettes as 
above and Taylor expand in powers of a. To be specific we note that the 1x1 
plaquette is the sum: 

Ax{x, y - 1) + ^y{^ + ^) ~ ^ + 1) ~ ^y^^ ~ f ' ^^'^^ 
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while the 2x1 and 1x2 plaquettes are the sums: 

A^{x, y-a)+ Ay{x + |, ?/ - |) + Ay(x + ^,y + |) 

- A^(x, y + a)- Ay(x - |, |/ + |) - Ay{x - |, y - |), (18) 

and 

Ax{x -^:y-^)+ A^{x + -,?/--)+ Ay{x + a, y) 

- ^x(x - |, y + |) - A^{x + |, y + |) - Ay{x - a, y) (19) 
respectively. The result of the expansion in a is: 

a ^Pii = (a^Aj, - a^/A^;) + — (a^^^Ay - dyyyA^) + O(a^) 

The results for the two remaining plaquettes are: 

a? 

a-^Pl2 = 2{d^Ay - dyA,) + —{d^:CxAy + Sd^yyAy - AdyyyA.,) + O {a') 

^2 

a ^Pai = '^{.dxAy - dyA^) + — (-^j^j^yA^. - 'id^xyAx + 4a3;^^Aj^) + O(a^) 
It is straightforward to verify the following: 

-2 2 2 «^ 2 2 4 

a = F + Y^(~(^a;a;^j/) ~ i^yy^x) + ^xyAxdxxAy + dxyAydyyAx) + 0(a ) 

«-'(A'2 + ^li) = - —[{dxxAyf + (a,,^.)^] + O(a^) 

g^2 

+ ~^[(c^xa;^2/)(c^a:2/^a;) + (^w^x) (^aiy^j/)] 

-a2[(9,,A,)2 + (6>,,^)^] (20) 

Now we may consider the linear combination found above and see what is the 
outcome: 

with the Taylor expansion: 

F' + Y2''^dxFf + {dyFf] (21) 

+ '^^[{dxxAy){dxyAx) + {dyyAx){dxyAy) - {dxxAyf - [dyyAxf] + O(a^) (22) 

We observe that several o? terms do not vanish with this (or any other) choice 
of parameters. We have a difficulty, stemming from the nature of the actual lattice 
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expression of the gauge fields. A possibility to eliminate these unwanted terms 
might be to employ further different kinds of plaquettes; however, what really 
happens is that, when bigger plaquettes are considered, new terms appear that 
cannot vanish against existing terms. 

An interesting remark is that one may choose different coefficients from the 
ones based on and find nicer expressions on the right hand side: 



We have decided to use the standard non-compact Wilson action (with no 
improvement) for the gauge part, mainly because of this difficulty and the con- 
sideration of the fact that we intend to use the action in the weak gauge coupling 
regime, so the subleading terms are not expected to be too serious. 

2.2 The gauge-scalar part 

We now go ahead with the gauge-scalar sector of the action. The part that needs 
improvement is of course the kinetic term, the only term involving derivatives. We 
feel that it is important to improve this part mainly, since its big autocorrelation 
times make mandatory a quicker approach to the continuum limit. Following the 
scheme of the previous subsection, we first consider the lattice embedded in a 
continuous space time and afterwards we turn to the actual problem that we face 
on the lattice. 

2.2.1 Lattice embedded in the continuum 

We start by writing down the continuum kinetic term in the fi direction for the 
scalar field: 

(j)*{d^ - iA^y(j) + h.c. = - (j)*(j)A^ - t(j)*(j)A' - 2i(j)*'(f)A + h.c, (23) 

where the primes denote differentiations with respect to and by A we under- 
stand Afj,. 

The kinetic term in the continuum involves the expression: 

0*(x)Per''"''^'"'">(x + a/i) + h.c. (24) 
Thus we choose to approximate this kinetic term by the expression: 

Shi = 0*(x)e''^/o '^*^("+"*^)0(x + ajl) + h.c. (25) 
and Taylor expand it in powers of a. The result is found to be: 

(f)*"(f) - (j)*(j)A^ - i(j)*(j)A' - 2i(f)* 
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+a\-^(j)*"(j)A^ - <P*'<pAA + + O(a^) + h.c. 

As we would like to eliminate the subleading terms, we can add next-to-nearest 
neighbour terms with suitable coefficients. We consider only up to second neigh- 
bours, that is we consider terms of the form: 

Sh2 = 0*(x)e'"/o + 2afi) + h.c. (26) 

The result of the Taylor expansion contains terms qualitatively similar to the 
previous ones: 

4(0*'> - 0*0^2 - i(j)*(t)A! - 2z0*VA) 
+lQa\-^-<P*<PAA' - ^<p*<pA' - \<P*<PA' + '-<P*<PA^A) 

+ 16a2(-^0*"0A2 - (jf'(j)AA + + O(a^) + h.c. 

It is easily seen that it is possible to choose the coefficients such that the o? 
subleading terms vanish. One need only consider the combination +^Shi — i^5'/i2- 
It is trivial to check that the result for the Taylor expansion of this combination 
reads: 

- 0*0^2 - i(j)*(t)A - 2#*>A + O(a^) + h.c, 
which is, actually the continuum action (pSj). 



2.2.2 Actual lattice formulation 



As in the previous case, on the lattice we don't have exactly the forms (^ , |26|) 
for the kinetic terms. The scalar field kinetic term before improvement reads: 
Sx./i 4>* {^)Uxfi4>{^ + O'f^) + h.c. To begin with, we write down the Taylor expansion 
of the expression 5*^"** = (/)*(a;)e*"^''^^"'"t/^V(3^ + c^A) + h.c. The result is: 

(f)*"(f) - i(f)*(l)A' - 2i(f)*'(f)A - (f)*(l)A'^ + h.c 

+a\-i^*(pA^A' - -(j)*"(j)A^ - (j)*'(j)AA') 
2 2 

+ a\-^^*(pA"A - U(P*'(j)A" - + O(a^) + h.c (27) 
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We see immediately that we have terms of order a?, which we would like to discard 
by the improved action. Our step towards the improvement will be to consider 
the next-to-nearest neighbor terms, namely J2x,fi(p*{^)UxfiUx+ii,ii(p{x + 2afi) + h.c. 
In the following lines we give the Taylor expansion of the expression 5^2** = 

0*(^^)giaA^{x+f A)giaA^.(x+f + 2a/i) + h.c. : 

A(f)*"(f) - Ai(f)*(j)A' - 8i(f)*'(f)A - A(j)*(j)A^ 

A 1 fi 

3 3 
W{Si<i)*(pA^A' - S<i)*"(pA^ - 160*VAv4') 

+a\-^-^i<j)*"'<j)A + U*""<P - Si<P*"<pA) 
o o 

+ a^i-h(l)*<pA'A - 5i^*'^A" - -i^*(PA"') + O(a^) + h.c. (28) 

6 

The next natural step to be taken is, of course, to use the coefficients found before 
and see whether the subleading contributions go away or not. The result for 
the linear combination +^Sl^l* — ]^5'^2** i^- 

(P*"<j) - i(P*(j)A' - 2i(p*'(j)A - (f)*(j)A^ + a^{^c^*(f)A"A+ ^i<p*'<pA" + ^i(f)*(f)A"') (29) 

We observe that the contributions do not vanish completely; the terms in the 
last lines of equations (pTj) and (^) survive. However, we observe that a good deal 
of the contributions (9 terms out of 12), present in the non- improved expression 
( |27D have disappeared. The conclusion is that we don't really manage to get rid 
of all the subleading contributions of order a^, but we expect that we approach 
closer to the continuum limit, so, presumably, the numerical behaviour of the 
improved action should be better; this has been confirmed by the simulations. 

Thus, gathering everything together, we conclude that the improved action 
reads: 

X 0<i<j X 0<i 

-^{^*{x)ip{x) - ip*{x)U^{x)Ui{x + iMx + 2i))] 
+ ^[(1 - 2Pn - 3lp,)^*ixMx) + PRiip*{xMx)f] (30) 

X ^ 

where Fij = A{Aj{x) - AjAi(x). 

The lattice parameters and the (three-dimensional) continuum ones are re- 
lated as follows: 

= 4 '''' 
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R 



4/3, 



20, 



2 1 — ^jPh 



2(3_ 



R 



y - {1 + Ax) 
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where S = 3.176 and S' = 2.752. In the appendix we prove the relation (|33|) 



(32) 
(33) 



3 Results 

We used the Metropohs algorithm for the updating of both the gauge and the 
Higgs field. It is known that the scalar fields have much longer autocorrelation 
times than the gauge fields. Thus, special care must be taken to increase the 
efficiency of the updating for the Higgs field. We made the following additions to 
the Metropolis updating procedure @]: 

a) Global radial update: We update the radial part of the Higgs field 
by multiplying it by the same factor at all sites: R{x) — >■ e^R{x), where ^ G 
[—e, e] is randomly chosen. The quantity e is adjusted such that the acceptance 
rate is kept between 0.6 and 0.7. The probability for the updating is -P(0 = 
min{l, exp(2V^(^ — A5'(^))} where AS'(^) is the change in action, while the 2VC, 
term comes from the change in the measure. 

b) Higgs field overrelaxation: We write the Higgs potential at x in the 
form: 

Viifix)) = -a ■ F + R\x) + pRiR^x) - if (34) 



where 



R{x) cosx{x) 
R{x) sinx(x) 



F = 



F2 



4 . . 1 

Fi = ^[-i?(f+i) cos{x{x+i)+e{x))-—R{x+2i) cos{x{x+2i)+e{x)+e{x+i))], 

I 

4 . . 1 . 

F2 - sin(x(x+i)+^(f ))-— i?(f+2i) sin(x(x+2i)+0(f)+0(f+i))]. 

I 

We can perform the change of variables: (a, F) — > (X, F, Y) , where 

. F 



:, X = a-f, Y = a-Xf. (35) 

The potential may be rewritten in terms of the new variables: 
V{X, F, Y) = -XF + {1 + 2(3r{Y'-1))X' + Y\1- 2(3n) + PrIX' + Y') . (36) 
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The updating of Y is done simply by the reflection: 

Y ^ Y' = -Y. (37) 

The updating of X is performed by solving the equation: 

V{X\F,Y) = V{X,F,Y) (38) 

with respect to X' . Noting that X' = X is obviously a solution, we may factor 
out the quantity X' — X and reduce the quartic equation into a cubic one, which 
may be solved. The change X — > X' is accepted with probability: -P(X') = 
min{Po, 1}, where P, ^ '-^^ r%'l'^r 

For our Monte-Carlo simulations we used cubic lattices with volumes V = 
8^, 12^, 16^. For each volume we performed 30000-50000 thermalization sweeps 
and 60000-100000 measurements. We have set the value of x equal to 0.0463. 
According to the relation (10) using mw = 80.6GeV and g = this value of x 
corresponds to Higgs field mass mn = SOGeV. For each value of Ph we determine 
the value of by using relation (|32D. The phase transition is expected to be of 
first order for such a low mass of the scalar field. 

We have used five quantities to determine the phase transition points: 

1. The distribution N{Eunk) of Bunk- 

2. The susceptibility of Eunk = Y.x,i ^* {x)Ui{x)Q{x +i) (we have set ip{x) = 
i?(x)e^^(^) = R{x)n{x) ): 

S{Eiink) = {ElinkY > — < Eiink >^)- 



3. The susceptibility of R2 = yJ2xR'^{^) '■ 

S{R2) = V{< {R2f > - <R2 >^). 



4. The Binder cumulant of Ennk- 

< {EiinkY > 



C{Eiink) — 1 



5. The Binder cumulant of i? = R{x): 

^^^^ - ^ ~ 3<{Ry 



The pseudocritical /9^(^, V) values have been found by determining (a) equal 
heights of the two peaks of the distribution N{Eiink), (b) the maxima of the 
quantities S^Eunk), S{R2) and (c) the minima of the cumulants C (Eunk) , C (R) . 
The values P^{A,V) depend on the specific quantity (denoted by A) which has 
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Figure 1: Susceptibility of Eunk 



been employed, as well as on the volume V. It has to be noticed that while 
searching we have made use of the Ferrenberg-Swendsen reweighting technique 
(11121) to find the pseudocritical Ph for the volume 16^. 

In figure 1 we depict the behaviour of the susceptibility Si^Eunk) versus j3h for 
three lattice volumes. The curves are fitted through the data; for simplicity we 
give the actual measurements for the largest volume. The curves represent the 
data quite nicely. In calculating the error bars we first found the integrated auto- 
correlation times Tint{A) for the relevant quantities A and constructed samples of 
data separated by a number of steps greater than Tint{A). Then the errors have 
been calculated by the Jacknife method ([0), using the samples constructed 
according to the procedure just described. Notice that the peak values increase 
almost linearly with the volume which is characteristic of a first order transition. 

In figures 2 and 3 we depict the behaviour of the Binder cumulants C{Eiink) 
and C{R) for three lattice sizes. Again, we show the real measurements for 
the largest volume. The error bars have been calculated, also, by the Jacknife 
method. The volume dependence of the cumulants display evidence for a first 
order phase transition. 

The different values of Pl{A, V) corresponding to the quantities A are due to 
the finite lattices used. So, we should extrapolate these values to infinite volume 
adopting the ansatz: 

V) = /3r(oo) + 

The extrapolated value Pf^^{oo) should not depend on quantity A because this 
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Figure 2: Cumulant of Eunk 
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Figure 3: Cumulant of R 
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Figure 4: Extrapolation for Pg = S 



is the infinite volume extrapolation for the critical point. 

Figure 4 shows the extrapolation to infinite volume using data for the pseudo- 
critical V) values obtained from the various quantities A versus the inverse 
lattice volume, along with the linear fits to the data. The error bars in (31{A, V) 
have been calculated from the statistical error of the values of the quantities A 
at the critical point. One can observe that, at, finite volumes, the smallest pseu- 
docritical values are given by the cumulant of Eunk', then, in ascending order, the 
values given by the cumulant of R, the equal height, the susceptibility of Eunk 
and the susceptibility of R2. Also, we notice that the infinite volume extrapo- 
lation is almost independent from the specific quantity used; the differences at 
the point = between the various extrapolated values are less than 10~^. The 
critical values lie in the interval (0.269176, 0.269183). In our previous publica- 
tion [|I| where we worked with the same model, but without any improvement 
to the Higgs part of the action, we had found that for the same value of the 
critical (3^ values were lying in the interval (0.336932, 0.336940). Although the 
precision is comparable, one should notice that the result presented here has been 
found by using three times smaller lattice volumes than the previous one which 
means much shorter computer time. Evidently, this result is due to the effect 
of improving the lattice Higgs action which provides a quicker approach to the 
thermodynamic limit. 

We can, now, predict the critical temperature Tcr- Actually, the quantity 
yields j/cr through equations ( P2| , P^); then equation (^) gives T^r = 130.64(19). 
In [U using 2-loop calculations it was found T^r = 131.18(14) while a 1-loop 
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calculation would give T^r = 130.74(14). Thus, the two 1-loop results are almost 
identical, but the first one has been achieved in a more economical way due to 
the improved action used. 
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4 Appendix 



We want to prove the relation (p3D, which relates the masses on the lattice and 
in the continuum up to 1-loop perturbation theory. In general we have: 

m| = ml — 5m\{h) (39) 

The counterterm will be calculated by considering the - terms from the 1-loop 
lattice effective potential ([|3)- 



We consider the pure gauge and gauge fixing part of the action (|30D . We use 
the relation Bn = and write down the kinetic term as: 

'-'9 g^a 




Sg 

Let us comment here a little bit about the gauge fixing term, which is the new el- 
ement here. For the perturbative treatment of Higgs models it is usual to employ 
the gauges. This choice is dictated by its simplicity, since then the mixing 
terms between the gauge field and the would-be Goldstone bosons vanish. How- 
ever, one should be careful with this gauge, since the so-called Nielsen identities 
[ P^ should be satisfied JT^ . This severely restricts the possible gauges; viola- 
tion of these identities will lead to unphysical results. This is the reason that we 
have chosen to stick to the Lorentz gauge fixing for the gauge action: 

Lgf = Y^T.T.[M^)-M^-^)r- 

This gauge does not face the previous complications, so we prefer to stay on the 
safe side, at the price, of course, of having also to deal with the non-vanishing 
mixed terms. Note that we take the Landau gauge ^ = at the end of the 
calculations. 

Going to the Fourier space (notice that A^{x) = J -^^e^^^^^^^ Af^{p) ), we 
find 



+ = ^ / ^ YIUpW' + ^,{\ - mP^U-p) (40) 

where 

2 . Pitt 

Pi = -sm— (41) 

Next we consider the part of the action involving the scalar fields. We follow 
essentially the same procedure as above, but a step that should be taken is the 
decomposition = :^(0o + 0i + i4>2) of the scalar field. Then the rescaling 

(fa)^*^ = has to be performed to get the corresponding part of the action in 
continuum form. The Fourier transform of the relevant part of the action reads: 
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+\ J + + ml)(j)2{-p) 

+ / |^^<W(^«*^(-P) + 1 1 ^s^MCn^r + ^PI)M-P) (42) 
with = jpi — jKpi. Let us collect some additional notations that have just been 



4^1' 12 

used: 



Pi = -szn— -, p =}_^p-, p =}_^Pi, (43) 



TT^T = 9l(l>l, rnj = mliiJ,) + 3X^(1)1, ml = ml{n) + Ag^Q. (44) 

Now, considering the quadratic part of the action, which is contained in the 
above equations, it is easy to read out the propagators Di, D2 of the fields 4>i, 02 
respectively: 

D^' = p' + ^p' + ml, D^' =p' + ^p' + ml (45) 

In the equation above the terms proportional to arise directly from the im- 
provement terms in the scalar sector concerning the next-to-nearest-neighbour 
contribution. 

The effective potential at one loop is found using the relation: 

Vl-^""" ^ -ln{Z), 

where Z = J[D(pi\[D(l)2][DA]e~^ . Keeping only the quadratic part of the action 
we have just Gaussian integrations, so we easily get the result: 

\ I (Cf^^^^^'^^"^^? " ^'^^^ ^^^^ 



- A^.- = if + ml + ^P.P.)5.- + ;^(r^ - l)m -'-Jf^, (47) 

iZ CL Qj -L/o 



where 

We write the gauge propagator in this form to display the contribution = 
^|r^, which is due to the mixing term between Ai and the imaginary part of 
the Higgs field. 

An integral which appears very frequently and should be calculated is the 
following: 

1 f^/^ d^p , ^ ^2 ~4 



If we differentiate the above with respect of m we take: 
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dl{m) = mK{m)dm (49) 

where 

with: = {amy, = AJ^iSin^^ and = 16J2iSin^Y ^^^^ point we 
follow (||15|) where the expansion of (^) in powers of M = am is to be used. 
In the following we denote by B the "Brillouin zone" [— 7r,+7r]^, in addition, 
TP = p"^ + j2p'^- Then the following equality holds: 

d^p 1 f d^p 1 f d^p 



B (27r)3 p2 + _Lp4 + J^4-2 (277)3 (n2)2 J b {2tt)^ {U^ + M^) 

The first integral equals the contribution of the second integral can be shown 
to be of 0{a), so it will not contribute to the infinite part. Gathering everything 
together yields: 

E' 

Kim) = h finite 

Ana 

Hence, integrating (^) we get: 

1 S' 

I(m) = + finite (51) 

2 Ana 



where 



is being calculated numerically at the value of S' = 2.752. 

Up to this point we are ready to calculate the infinite part in the effective 
potential (behaving like -) that is due to the scalar fields only. The mass terms 
coming from the fields 0i, 02 are: 

IE' 2 IE' 2 

m,, m, 

2 Ana ^ 2 Ana ^ 

respectively. Recalling (0), their derivatives with respect to the classical field 
are: 



a2 



The next step is to calculate the gauge contribution to the effective potential, 
that is the third term in equation (^6|). The calculations are straightforward, 
but quite tedious; we just mention that we directly expand the determinant and 
keep contributions only up to order a^, since the terms will yield merely finite 
results, which are not our concern here. We only write down the final result: 

^2"'°"' = / ^sMP' + mlf + (/V - inO 
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Plpl+PlPl+PlPl. , ^/ 4n 



+ Inil + 2 ,^ ' ' "^"^ ) + 0(a^)] (54) 

We note that the second term is exactly the same as the one appearing in the 
continuum counterpart of the model, so they cancel upon comparison with the 
continuum model. These terms do not depend on the mass, so they present no 
interest for the calculation of the effective potential anyway. We should note that 
we have not given the full result of the third term; to keep things simple, we gave 
its value for = 0. The above expression can be written as: 

The integrals appearing in equation (|55| ) have already been computed; the infinite 
part of the effective potential reads: 

-mj. + - — my + —r-mrp) (56) 
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We invoke equation (|4^ ) and calculate again the second derivative with respect 
to the classical field 0o to find out that the infinite part reads: 

^3/ S 1 N 

T^i;^ + 4^^+12^^ ^^^^ 

Collecting all the infinite contributions we can match the counterterm equation 
and this yields: 



^ Pft 4 47r 47r 12 



2/5.(75- -3T-^)=y-(l + 4x) 
Pft 4 

which is the equation (|33D . 
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